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Andreev reflection between a normal metal and a superconductor whose order parameter exhibits 
quantum phase fluctuations is examined. The approach chosen is non perturbative in the tunneling 
Hamiltonian, and enables to probe the whole range of voltage biases up to the gap amplitude. Results 
are illustrated using the one-dimensional Josephson junction array model previously introduced in 
the linear response regime. Phase fluctuations are shown to affect the differential conductance and 
are compared to the result of Blonder, Tinkham and Klapwijk for a rigid BCS superconductor. The 
noise spectrum of the Andreev current is also obtained and its second derivative with respect to 
frequency is proposed as a direct tool to analyze the phase fluctuations. 



I. INTRODUCTION 

In the last decades, a considerable effort has been devoted towards the study of the transport properties of nor- 
mal metal-superconductor (NS) junctions. The situation for the Andreev current and finite frequency noise is well 
understood when the superconductor is of the BCS type. More recent works have dealt with superconductors whose 
order parameter has a d-wave symmetry [[l]|^] . The role of collective modes arising from the phase of the fluctuations 
of the superconductor has also been addressed in the framework of linear response theory in high-T c materials [gj, 
as well as in one-dimensional array of Josephson junctions Q. Some recent attempts have also included the effect of 
classical phase fluctuations of the order parameter on Andreev transport, either in the tunneling regime Q or using 
Bogolubov-de Gennes equations Q. 
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FIG. 1. Andreev reflection (bias voltage eV = fiN — Ms): a) In a rigid superconductor, the two electrons (one electron e~ 
and a time reversed hole h + ) have opposite energies ±hui with respect to the superconductor chemical potential us, and a 
Cooper pair is created on the right hand side, b) In the presence of phase fluctuations, the superconductor absorbs a collective 
mode with energy h£l, so that the energies of the two incoming electrons are no more symmetrical with respect to fj,s- 
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A physical picture of the Andreev process is drawn in Fig. 111. A voltage bias is applied, such that the chemical 
potential of the metal /in lies eV above that of the superconductor //g. The origin of the energies is taken at us- 
In Andreev reflection with conventional, rigid superconductors, Fig. [l]a, an electron below /ijvj with an energy hu> is 
reflected as a hole and a Cooper pair propagates in the superconductor. This pair has zero energy. The reflected hole 
can be interpreted as an electron below /is, with energy —huj, which enters the superconductor. At zero temperature, 
the only possible values of tko are between — eV and eV, see Fig. |l|a , because both electrons states must have an 
energy below /iAr. 

Turning to the case of a superconductor with phase fluctuations, Fig. |l|b: exciting a collective mode in the 
superconductor costs an energy Ufl. One can therefore take an electron of energy 7i(f2/2 + and reflect a hole, which 
is equivalent to a missing electron of energy— U{oj — £1/2). The allowed values of %ui must be smaller than eV — Ml/2 
(no electrons present above /xjv)- Note that here, one has to sum over all possible values of hil £ [0; 2eV], taking care 
of the fact that this is a weighted sum, with a probability distribution -P(il) which is determined by the Hamiltonian 
describing the phase. 

In this paper, the Andreev current is calculated non-perturbatively in the tunneling Hamiltonian for N-S junctions 
whose Hamiltonian includes quantum phase fluctuations. The approach is inspired from the Keldysh technique of 
references , and enables the calculation of the current- voltage characteristics of the Andreev current when the bias 
voltage eV is close to (but smaller than) the gap amplitude | A|. In particular, this corresponds to the experimental 
situation encountered in scanning tunneling microscope (STM) experiments performed on high-T c materials [fl0|-[l2[ . 
Although the calculation is restricted here to an s-wave order parameter, qualitative features can be derived for 
transport in N-S junctions for arbitrary superconductors, provided that there is a gap. Here, a phase-only effective 
Hamiltonian will be used, which was previously derived in the literature jL3|. In one dimension, this model is 
isomorphic to an array of Josephson junctions, with Josephson coupling energy Ej and charging energy Eq. The 
current and the finite frequency noise spectrum S(u>) are both derived, which requires to go beyond linear response 
theory. In particular, it will be shown that the second derivative of the noise with respect to the frequency d 2 S/dui 2 
gives a direct access to information about phase fluctuations, independently of the model chosen to describe the latter. 

There are other situations where Andreev scattering is mediated by excitations at the normal supraconductor 
boundary. If the normal side is replaced by a ferromagnet ||, the spin of electrons pairs emerging from/entering 
the superconductor can suffer a spin flip accompanied by the emission/destruction of a magnon. This additional 
scattering channel effectively may enhance the Andreev current. Except for the case of a half-metal, electrons with 
opposite spins as well as electrons with the same spin then contribute to the Andreev current. The differences with 
our phase fluctuations calculation will addressed below. 

The paper is organized as follows: in section [il] , the model Hamiltonian is discussed. In section [II, the perturbative 
scheme which allows to derive the relevant Green's functions for the transport properties is described. In section IV, 
the phase fluctuations are introduced in this framework. The current- voltage characteristics associated with the linear 
Josephson junction array model of Ref. Q for the phase appears in section [v|, as a function of the physical parameters 
Ej, Eq, as well as the barrier transparency. 



II. MODEL HAMILTONIAN 



The Hamiltonian describing the NS junction is composed of three terms: H — Hn + H t + Hjja- The first term is 
simply the Hamiltonian for the metal which is specified to a one band model Hn — ^ k efc,crcl a Ck t „ with width W. 
The second term H t is the tunneling term, 

H t = rcjv, CT c| CT + h.c, (1) 

N,S,a 

where T is a hopping amplitude which transfers electrons to/from the superconductor. 

In the superconductor, it is assumed that the low energy lying excitations are given by fluctuations of the phase of 
the order parameter. If we forget the poorly screened long-range part of the Coulomb interaction, the Hamiltonian is 
the same as the one describing an array of Josephson junctions p|,|l4|]. 

Hj JA = ]T Ej COS& - Bj) + £ ^ , (2) 

where Ej is the Josephson coupling (between islands i and j), the operator 8i is the phase of the order parameter 

2 _ 

at "site" i, qi is the charge on island i and ^ is the charging energy of island i. Eq. (0) can be put on a firmer 
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basis, starting from microscopic models. In particular, a phase-only Hamiltonian H$ has been derived from various 
microscopic models such as the attractive Hubbard model p3[ : 

h 2 Na ,„ „ . 1 m*v 2 9 x -> 7re 2 p;Oj 

^ = E-2^r «»(** - «w + E 2^ + E — i^TT ' (3) 

i,j i 2)J 

where AT§ is the bare superfluid density, e is the dielectric constant, | r< — fj | is the distance between site i and j. 
The first term is the Josephson coupling between location i and location j. The second term is the local density 
fluctuation energy (p — (p)) 2 /2xo, where (p) is the average density and %o a density susceptibility. The last term is 
due to poorly screened long range Coulomb interactions. Since electron and hole-quasiparticles are not quasiparticles 
for the superconductor, the injected electrons and holes will decay and form pairs that will excite the eigenmodes of 
H s . 

III. NON-EQUILIBRIUM GREEN'S FUNCTIONS AND TRANSPORT 

A voltage bias is introduced between the normal metal and the superconductor and the average current reads 

I(t) = jJ2 T ( c NA t ) c sAt))~h.c. (4) 

a 

Tunneling occurs from the bulk of the normal metal to the bulk of the superconductor through a narrow constriction, 
the tip of a STM for example. 

The Keldysh formalism is used to express the current in terms of Green's functions. The definitions and properties 
of this system have been studied in detail in Ref. Q . Here we adopt the same notations and start from their expression: 

2e f°° 

('(*)> = 72 l r | 2 / [9^Nii(t,h)Gs+ n (t u t) - g^+^^G+s^t^dh. (5) 
n J-oo 

The dressed Keldysh Green matrices on the superconducting side are expanded using the Dyson expansion 

G^(t,t') = [cT ss t^tG% 8 ] (t,t'), (6) 

where E is the self energy matrix describing the hopping from the normal side to the superconductor, and £i describes 
the opposite process. For simplicity, integrals over time have been omitted. Note that only subgap tunneling is 
considered here, allowing to set g^ = 0. Inserting the Dyson expansion, the current reads: 



(7) 



Op r°o roo rco 

<!(*)> = 74 llf / / / G r ssl2 (h,h)G a SS21 (t 3 ,t)x 

The dressed, retarded Green's function on the superconducting side is given by: 

G r S si2(t',t)=g r ssl2 (t',t)+ T (gl s {t'M)^ r s(tiM)G r ss {t 2 ,t)) dt 1 dt 2 + ..., (8) 

with the dressed self energy Eg defined as 

£ r s(h,t 2 )= I I t^tut'^x&AW'iMW^ , (9) 



OO pOO 



and the unperturbed momentum averaged Green's function in the superconductor reads : 

- wTPRrar"" s 5W " )e * (10) 

where a is an infinitcsimally small positive quantity and 8q the phase of the superconductor, which appears also in 
the gap parameter: A = |A| e l6 °. 
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IV. TRANSPORT IN THE PRESENCE OF PHASE FLUCTUATIONS 



In the presence of phase fluctuations, the Green's function on the superconducting side contains two time scales 
|^5| . First, fast time scales (of the order Ti/A) which were described in the previous section in Eq. jl0|). Secondly, 
"slow" time scales associated with the fluctuating phase 9(x, t): 

A = | A | e i<?(=M). (11) 

The phase of the order parameter has both spatial and temporal variations and is no longer equal to 6q = #(0,0). 
Tunneling is supposed to occur close to the interface. The typical dimension parallel to the interface will be small 
with respect to both the classical correlation length Jl6| and to the amplitude correlation length £ a . In high-T c 
superconductors, and £ a are of the order of a few nanometers. £5 is believed to be generally not much larger than 
£ a , and £ a is much smaller than in BCS superconductors. We shall thus replace x by zero. The temporal variations of 
are due to quantum fluctuations; the phase of the order parameter and the number of pairs are conjugated variables. 
These variations occur on time scales of the order at least of the plasma frequency io v = \Z8EqEj /h, with Eq = <7q/2C 
and Ej appear in the definition of the Hamiltonian of Eq. (^|). Equivalently, huj p <| A |, which is the contrary to 
what usually occurs in BCS superconductors. The dynamic nature of the phase is reflected in the anomalous Green's 
functions fx: 

5s'5i2^ ) t)=PS5i 2 (*'^)e ii ^ (0 ' t)+e{0 ' i ' ) ) (12) 

where 5<j'<j 12 (i', t) is the inverse Fourier transform of 95512(^)1 an d 9ss2i(t' > nas ^ ne opposite phase. Because of the 
slow temporal variation of the phase, the exponential in Eq. (fl2|): exp[i(0(O, t) + 0(O,t'))/2] « exp[i6(0, t')] without 
affecting the results in a significant manner. 



A. Current and differential conductance 

Substituting Eq. ( |l2| ) in the current, the fourth (lowest non vanishing) order term in T which contributes to the 
current is: 

c\ fOO />oo />oo 

(Ht)) = ^\r\ 4 / / 3552i(i3,t)5S5i 2 (^i^2)(e^)e-^)) 

h J-00J-00J-00 (13) 

x \9nnu (*' t i)9NN22^h- t h) - 9NNn(t, *i)ffjv-jv22(*2) tafjdtidtzdta. 

The bracket (...) signifies that an average has been performed over the dynamical degrees of freedom of the phase. 
Depending on the nature of the phase, this will correspond to a Gaussian model - the so called ordered phase - or to 
a nonlinear model |L7| (see below). Note that this expression is analogous to the one derived several decades ago by 
Combescot et al. for inelastic tunneling ]18[ |: there the working assumption was that the inelastic coupling did not 
couple directly the two electrodes. Here, this same assumption is natural because the phase fluctuations happen only 
on the superconducting side. The correlator of the phase and its Fourier transform are denoted: 

p(t) = ( e <9«M) e -«(0,0)} ) (14a) 
1 f°° 

P(u) = — / e^p(t)dt . (14b) 

P(uj) can also be viewed as the occupation probability of the collective phase modes with energy Tllu. Using the 
Fourier representation, the current becomes: 

r 00 poo 



(i(t)) = -± \r\* / flS OT1 (w + n/2)§s sl2 (o;-n/2) 

fill J— QOJ — OO 



- 00 «/ - — oc (15) 

x 9nnii(^ + n / 2 )9 N 1v22(u - n / 2 ) - gNNiii^ + n / 2 )9NN22( UJ - n / 2 ) P(tydnduj. 

Next, for frequencies below the gap, one can approximate gs'siji^ ^ ^/2) — g'ssij^ LjJ ) — ^/ ^H^^/l — h 2 uj 2 /\ A| 2 ^ (W 
is the bandwidth), and the current becomes: 
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<^)> = 7-3 |r| 4 / / ~g a ssM~9 r ssM , s 

nn J-ooJ-oo (16) 



where the two causal Green function elements in Eq. ( [L6| ) are related to the electron density of states p(tuv): 

StvniiM = 27rf/3n (hu - eV)n F (huj - eV) , (17a) 
9nn22( uj ) = 2nip22(hu + eV) [n F (huj + eV) - l] , (17b) 

with pu(huj — eV) = p22{— nui + eV) = p(hu) ~ H/ttW. 

Using the property g%S2ii L0 ) — [ffssi^C^)]* > the current in the absence of phase fluctuation is recovered by setting 
P(fi) to be a delta function: 

o I "p |4 /"Oo 

{I)o = JJi W^J I #?S2iM I 2 [nj?(fiw - el/) - n F (huj + eV)\dw, (18) 

where subscript means that in this limit Eq/Ej = 0. The known result for the Andreev current is recovered [pj. 
It is well known however that even in the absence of phase fluctuations, a perturbative calculation to fourth order 
in r is not sufficient to give a satisfactory answer for the whole range of voltage biases between and A. Indeed, a 
resummation procedure, explicited by Eq. (^|), has to be carried out. In the absence of phase fluctuations, it is thus 
sufficient to replace gss2i( U! ) by G r SS21 (uj): this allows to recover the scattering theory results of BTK JTgj ] . 
Turning back to the lowest order contribution of the current in the presence of phase fluctuations, one obtains: 

o I p 14 poo i-oo 



hhw* J J x — — " (19) 

x n F (huj-eV + hD,/2)-n F (huj + eV -HQ/2) P(n)dflduj. 

Note that -P(fi) decreases when f2 becomes large. In particular, the integrated contribution J„ P(Q')d£l' becomes 
negligible well before hQ > 2eV, allowing modify the upper bound of the integral over tt. Simultaneously, the advanced 
(retarded) Green's functions are dressed as before and take the form : 

G ssi 2 M = p^rx fssM' (2°) 
so that at zero temperature the current reads: 

16elrl 4 f 2 * /-Tf-f i 



with 



DU) [\^-^D r (uj)D^)=(^-^)\Lr^f 

lAi 2 ,\ / in 4 \ 2 .in 4 



r? U )[ 1 + W) +4 W^ (22) 

In a more rigorous approach, the expansions of G r ssl2 {to) have to be performed in the presence of phase fluctuations. 
This means that they should contain not only linear terms in P(f2) but also higher order correlators of the exponenti- 
ated phases. Here, in the expression for the current in Eq. (pl|), only the first order in P(f2) has been retained, which 
constitutes the analog of a weak inelastic coupling - single phonon - approximation of Ref. p8| . 



The differential conductance d(I)/dV can be straightforwardly computed from Eq. (21) 

d(I) 16e 2 |r 



A l / P(y) „, eV v , dy. (23) 



dv hh 2 w± 1 1 Jo £>(ir-f) 
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B. Finite frequency noise 



Another useful tool which probes the fluctuations is the finite frequency noise. The symmetrized real time noise 
correlator is given by: 



S(t - t') = (51 (t) 51(f)) + (5I(t r ) 5I(t)) 



(24) 



where 5I(t) = I(t) — (I(t)). This current-current correlator is expressed in terms of the Keldysh Green's functions: 
using the self energy matrix, we can write it as a trace. 



S(t, 0) = -e 2 Tr [ { [a z ±G+ N (0, t)<7 z EGg + (t, 0) - a z ^G+ N (0, t)a z £6g+(t, 0)) + ( t «-> ; S «-> N ) } 



+ { (0,t)«-(t,0)} 



(25) 



where a z is a Pauli matrix. 

The lowest non vanishing contribution (order 4 in T) is extracted by expanding the Green's functions in Eq. (|2 
The corresponding Dyson equations read (time integrations are implicit): 

&SN = GsS^fiffN ' 

NN = 9nn i (26) 

where we have purposely written the Green's functions to the lowest order which they contribute. Inserting these 
expressions in the real time correlator gives: 

/+oo r+oo 
dh / dt 2 Ti { a z ±G ss (0,t 1 )±^g+ N (t 1 ,t)a z tG ss (t,t 2 )t^ N + N (t 2 ,0) 
-oo J — oo 

+a z ^g+ N (0,t 1 )tG a ss (t 1 ,t)a z ^g N +(t,t 2 )m% s (t 2 ,0) 
-a z £6 ss (0, h)±*g+„ (ti, t 2 )tG a ss (t 2 , t)a z ^g-+ (t, 0) 

-a z Et.9+ w (0, t)*,£(% s (f, ti)St^+ (ti, t 2 )±G a ss (t 2 , 0)} + {o ^ t}] . (27) 

The procedure for including the phase fluctuations in the noise correlator is identical to that used in the expression 
of the current of Eq. (R). However, here additional contributions which are proportional to G r ss G r ss or to G SS G SS 
occur, but these turn out not contribute to the noise to lowest order in the phase correlator. 

The current noise spectrum S(u>) is the Fourier transform of S(t). In the subgap regime, the only non-vanishing 
elements of the causal Green's functions are the off diagonal ones : g^g^ — 9s's2i — h/W. The current noise spectrum 
can be expressed in terms of the off diagonal Keldysh Green's functions: 



S(u) 



4e 2 



dnp(n) 



du>i 



h 2 2irW 2 jo 

The integration over u)\ defines the relevant energy intervals which contribute to the noise: 



(28) 



16e 2 irl- 



2tt W i 



dnp(n)e 



2eV 



uj — Q 



/2eV 



LU — il 



2 ^+" (2eV 

dnP{Q) 1—+UJ-Q 



(29) 



where 8(f2) is the Heaviside function. 

In the limit of low frequencies, the noise becomes: 



G 



S(u = 0) 



_ 32e 2 |r| 4 
= 4e(7), 



(30) 



which corresponds to the Schottky formula. Note that this constitutes an illustration of Schottky formula for a 
situation where the charge transfer is effectively inelastic: as shown in Fig. [l]b, the transfer of two electrons in the 
superconductor generates phase quanta and thus constitutes an inelastic process. These processes are included when 
the distribution P(u>) in the current and noise deviates from a delta function. 

The current noise spectrum can be measured experimentally and provides a direct way to obtain information on 
P{uj) and thus on phase fluctuations. In the tunnel limit, the distribution P(lu) is related to the second derivative of 
the current noise spectrum by the formula 



d 2 S 
duj 2 



V 



p 



2eV 



P 



2eV 
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2eV 



(31) 



The prefactor stands for the current in the absence of phase fluctuations :(/)o = (16e 2 V r |r| 4 )/(?iVF 4 ). 

Note that so far, no specific model for the phase fluctuations has been specified, as the only assumption used what 
the fact that -P(Sl) falls sufficiently rapidly under the gap. In the next section, current voltage characteristics are 
plotted using a specific model to describe the phase dynamics. 



V. APPLICATION TO THE ID JOSEPHSON-JUNCTION ARRAY MODEL 



One spatial dimension is specified, motivated by the fact that tunneling typically occurs near the tip of a scanning 
tunneling microscope. The dynamics of the phase can then be modeled by a one-dimensional Josephson-junction 
Array (JJA), which was previously studied in linear response g. The phase correlator in Eq. (14a) is then known. 
The parameters of this model are the Josephson energy Ej and the charging energy Eq (which in turn define the 
plasma frequency), together with the parameter k which depends on the ratio of the junction capacitance divided by 
the capacitance of the islands. If the latter capacitance dominates over that of the junction, k = 2. 

The importance of phase fluctuations is monitored by the ratio Tiu> p / Ej. In particular, for Twjp/Ej < tt/2, the 
phase Hamiltonian can be mapped to a Luttinger liquid model and P(f2) can be derived in a standard way from 
harmonic oscillator correlators. On the other hand, for hui p /Ej > 7r/2, the phase correlator decays exponentially in 
time. This is identified as the disordered phase. At hui p /Ej = tt/2, there is a Kosterlitz-Thouless transition between 
the ordered and the disordered phase. 



A. Ordered phase 



At zero temperature, the time dependent phase correlator is given by: 

r*»>p\t\ l-cosx dx 
la 



J V2 

p{t) = exp<^ W — 

n \ Ej 



x a/1 - (x/nupt) 2 



i sgn(t) 



sin x 



dx 



x a/1 - (x/KUJ p t) 2 



(32) 



with sgn(i) the sign function. Its Fourier transform P(Cl) is identically zero for O < 0. For > 0, it has a power law 
behavior close to £1 = 0: 



p(n) 



(33) 



and goes rapidly to zero for > nui p . 

For f2 = KLu p , P(O) shows a weaker singularity JpJ and diverges as (uj p -Q)- a , with a = 1/2 - tt- 1 y/2E / Ej (here 
this corresponds to fixing the parameter of Ref. ETk = 2). 

Turning to the transport properties, the effect of the phase fluctuations on the noise are probed. For convenience, 
we compare our results to those derived with a BTK model ]l9| , which corresponds to an interface with a rigid 
superconductor. The quantity 7 = T/W corresponds to the transparency of the barrier. 7 enters the prefactor of the 
delta function potential of BTK theory with the dependence T% 2 kFi^ — 7 2 )/2m7 (&f is the Fermi wave vector). 7 = 1 
corresponds to a perfect contact, while 7^0 describes an opaque barrier. 
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FIG. 2. Second derivative of the symmetrized noise spectrum versus frequency (in units of 4(e 2 / h)h-y 4 /u> p ) for different 
values of hu p /Ej, hujp/Ej — 0.01 (solid line), hu p /Ej = 0.2 (dashed line), %u p /Ej = 1 (dashed dotted line). The plasma 
frequency uj p is set to \A\/2h. 

In Fig. H the second derivative of the symmetrized noise spectrum d 2 S/du> 2 is plotted as a function of u>, for a 
perfectly good contact. Several values of the parameter fru> p /Ej are considered, starting from weak fluctuations to 
substantial collective excitations within the ordered phase. To illustrate the effect of phase fluctuations, the choice 
huj p = |A|/2 is made, so that all the features in P(fl) are located within the superconducting gap. 

In the absence of phase fluctuations, the finite frequency noise has a singular derivative at fiuo = 2eV |2(J. In Fig. 
|2|, this would imply a delta function peak in d 2 S/du> 2 . For low voltages, there is no deviation with respect to the 
results of scattering theory. The effect of phase fluctuations is twofold: first, the delta function broadens to a power 
law divergence (partially cut on the figure), and second, it acquires a secondary peak - previously discussed for P(fi) 
- at 2eV — nhujp. This illustrates the result of Eq. (|3l|): an experimental measurement of the current noise spectrum 
constitutes a direct diagnosis of the importance of phase fluctuations. 

Naturally, the differential conductance d(I)/dV in Eq. ( p3| ) is also affected by the phase fluctuations. It is studied 
below for different regimes: first, it is plotted as a function of bias, for both the case of a perfect contact and for the 
case of a weak transmission; second, the transparency is varied while the bias is fixed to a large value in order to 
explore the deviations from the tunnel limit. In each case, curves are obtained for several values of the ratio hoj p /Ej 
which characterizes the importance of the fluctuations. 

The differential conductance of a perfect junction is displayed in Fig. || as a horizontal line at d(I)/dV = 4e 2 /h 
In the presence of fluctuations, it saturates to the BTK value for biases larger than Tiu p /2 (here the bias is 

varied from to |A|): for sufficiently large eV, Q &V ^ P(£l)dfl is essentially equal to unity because for 7=1, D(x) 



which enters Eq. (|23j) is a constant. Phase fluctuations have a tendency to decrease the differential conductance at 
low bias. In particular for Tiuj p / Ej > 1, the significant contributions from the integral J^ eV ^ h P(£l)dSl are the cause 
for the deviations from the ideal conductance at large bias. 



4. 
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FIG. 3. Differential conductance d{I)/dV in units of e 2 /h versus the ratio eV/Tiu> p . A perfect contact (7 = 1) is considered 
and different values of hui p /Ej are displayed : hu) v / Ej = (BTK curve, solid line), hui p /Ej = 0.2 (dashed line) and Tiuj p /Ej — 1 
(dash-dotted line). 

The case of 7 = 1 is somewhat academic because perfect contacts are difficult to achieve in practice. The same 
curves are plotted now for 7 = 0.4 in Fig. ||. For intermediate biases, phase fluctuations do not affect the differential 
conductance drastically: at eV = hui p for instance, d(I)/dV is essentially constant (of the order 4e 2 j 2 /h), and all 
curves can not be distinguished. However, deviations occur both at small and large voltages. For small voltages, 
according to Eq. (|^) and due to the fact that P(0) = the differential conductance in the presence of any phase 
fluctuation is required to vanish (this is not really obvious in Fig. |I] because of the choice of scale). More dramatic is 
the fact that phase fluctuations play a role at large voltages, contrary to the high transparency regime. 




1 ■ eV/ficfl 



P 

FIG. 4. Same as Fig. for a transparency 7 = 0.4. 

In Fig. ||, the transparent regime is examined with a bias voltage eV — |A|. In contrast to the previous curves, 
the differential conductance is plotted versus the transparency of the barrier. The plasma frequency and the ratio 
hjJ v jEj are chosen as previously. 
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0. 0.2 0.4 0.6 0.8 1. 
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FIG. 5. Differential conductance (in units of e 2 /h) as a function of the transparency, for a voltage bias eV = |A| (same 
convention for hco p /Ej as in Fig. 3) 

For barrier transparencies smaller than 0.1, the differential conductance has power law dependence {d(I) / dV ~ 7 4 ). 
As the transparency is increased, in the presence of phase fluctuations an inflexion point appears and the rise of 
d(I)/dV is then slower. Note once again that it is reduced from the BTK case when phase fluctuations become 
important. As indicated previously in Fig. the differential conductance saturates to the BTK value at high 7 
regardless of the degree of phase fluctuations. 
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B. Disordered phase 



In the disordered phase (Uu> p /Ej > ir/2), the behavior of the Fourier transform of the time dependent phase 
correlator close to the transition has been previously derived from the XY model |17|. For small f2, 



P(fi) ~ e(n - UpC 1 )^ - UpC 1 )^ 1 ^, 



(34) 



where &(x) is the Heaviside function. The results of linear response theory of Ref. Q are briefly recalled and 
compared to the predictions of the Keldysh calculation of Eq. (pl|). In this phase, a threshold voltage Tiuip^ 1 has to 
be reached in order to obtain a non zero current. Here, the Kosterlitz-Thouless correlation length (in dimensionlcss 

b I 2 /it — \J Ej /8Eq ) , where b is a positive constant. Because the goal is to predict the 



units) reads £ = exp 

current as a function of voltage for both the tunnel limit and the transparent regime (for < eV <| A |), P(Q) 
needs to be characterized. P(Q) can be set to zero for Q larger than klo p . This assumption is reasonable below the 
Kosterlitz-Thouless transition. In accordance with Eq. (J34l), we adopt the form 



Using Eq. (^l| ) and (35), calculations can be performed analytically. 

For eV < hujp^ 1 /2, the average current vanishes. For < 2eV/hu>p < k, the average current reads : 



(35) 



(I) 



I67 4 



hh (1 + 7 4 ) 2 



(i-«r x ) 



nujp ^ 



3/2 



(36) 



where 



A = 




1 arctan 




(37a) 



(37b) 



where A is a scaling variable which varies between (perfect contact 7 = 1) and 1 (7 = and a plasma frequency 
ujp = |A|/?i). Note that the function / does not vary substantially between these two transparencies: f(A = 0) = 2/3 
and / (A = 1) ~ 0.71. The dependence of the current on the parameters 7, eV, aj p and £ is given by the prefactor in 
Eq. A 



For nhujp/2 < eV <\ A|, 



(I) 



where 



I67 4 



M(l+7 4 ) 2 2 



(38) 




- 2 



y(2-x) 



arctan 



For 7=1, the behavior of (i) is essentially the same as that of linear theory 
plotted in Fig. (with the choice uj p — A/4?i). 



(39) 



y(2-ar) y y 
for arbitrary w p . These results are 
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o. 0.5 1. e v/KRco 2 - 

p 

FIG. 6. Current (7) versus bias voltage in the disordered phase for hu> p /Ej — 2ty [in units of (4eaj p /7r)7 4 /(l+7 4 ) 2 ]: perfectly 
good contact (solid line) and for 7 = 0.4 (dashed line). 

For eV larger than Khui p /2, the curve (I) versus eV is a straight line. For less transparent barriers (7 ^ 1), two 
behaviors have to be distinguished. First, if the plasma frequency is smaller than the gap parameter, then the scaling 
variable A in Eq. ( |37b| ) is small and linear response is again retrieved. Second, if Tiu p is comparable to | A |, (I) 
increases typically as (2eV — nfiuip^ -1 /2) 3 / 2 for £ -1 /2 < 2eV/nhujp < 1. Then, for eV larger than Kfuv p /2, the 
response is non-linear. In particular, when the transparency is low, (I) increases faster than {2eV/ kT%u> p — £~ x /2) 3 / 2 
for eV > nhtOp/2. Overall, no sensible deviations from linear response are expected in the disordered limit. 



VI. CONCLUSION 



A formulation of quantum transport in normal metal superconductor junctions is developed in the case where the 
order parameter undergoes phase fluctuations, both in the tunnel limit and in the transparent regime. The Keldysh 
Green's function was used, allowing to go beyond linear perturbation theory, provided that the subgap regime is 
specified (this latter assumption enabled us to neglect some contributions in the Dyson expansion) . While we applied 
our results to the ID Josephson junction array model of Ref. Q, the result is quite general as any model for the 
phase correlator can be employed to calculate the current and the finite frequency noise. An interesting issue is that 
the zero frequency noise satisfies the Schottky formula (as should be), with a current describing inelastic processes. 
A direct correspondence between the phase correlator and the second derivative of the noise spectral density with 
respect to frequency has been established. For a highly transparent interface, phase fluctuations are shown to affect 
only the differential conductance at low voltages: at high voltage the I(V) characteristic is essentially that of a rigid 
superconductor, which can be described with the BTK model. On the other hand, for a less transparent barrier, 
deviations with respect to the BTK are shown to occur both at high and low voltages. This effect becomes more 
dramatic when the barrier is opaque. 

The present system bears similarities with other situations where Andreev reflection is accompanied by the emission 
of a bosonic excitation. As mentioned in the introduction, Ref. Q dealt with a ferromagnet whose magnons can flip 
the spin of one of two same spin electrons in order to contribute to the subgap conductance in a s-wave superconductor. 
Therefore, electrons with opposite spin contribute to the current as usual, but a new scattering channel - due to this 
one-magnon process - is opened, and enhances the Andreev current overall. The contrast with our calculation is that 
here, the current is always reduced by the excitation of phase modes. For a specified voltage bias, the more phase 
fluctuations, the more energetic are the modes which propagate in the Josephson junction array. On the side of the 
normal metal, this requires that the sum of the energies of the electron and hole has to be large. Yet the Fermi level 
restricts the impinging electron energy. This reduction of the available phase space for inelastic assisted Andreev 
reflection thus justifies a weaker current when the fluctuations are large. Also, note that the phase correlator which 
was used in this calculation typically appears in situations where a quantum mechanical degree of freedom is coupled 
to an electromagnetic environement pl| . Overall, if one includes multi-magnon processes to the situation of Ref. pj , 
both ferromagnetic/superconductor and normal- metal/ fluctuating superconductor junctions constitute illustrations 
of quantum dissipative systems ]22] ] , where the quantum process to be probed - Andreev reflection - is coupled to a 
physical environment (a bath of magnons/ "phase oscillators"). 
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Potential applications to physical systems are now addressed. While typical BCS superconductors have essentially 
amplitude fluctuations (except close to the superconducting transition) , the high-T c cuprates have a much lower phase 
stiffness and hence, are expected to show strong fluctuations of the phase of the order parameter p3| . Here, an s-wave 
order parameter was assumed for simplicity, whereas these materials are mainly d-wave, but our calculation could in 
principle be generalized to this latter situation. The two main differences between s-wave and d-wave materials are 
the following: first, in a d-wave material, there are nodes in the order parameter in specific directions, which results in 
a finite density of states for quasiparticles near the superconductor chemical potential. Second, low energy Andreev 
bound states may exist at this energy. Depending on the orientation of momentum, an electron-like quasiparticle and 
a hole-like quasiparticle in the superconductor see an order parameter with the opposite sign B. This generates a 
zero-bias conductance peak (ZBCP) which is a hallmark of a d-wave order parameter. 

At first sight, and for ideal interfaces, the connection of between our s-wave model and a d-wave compound could be 
made by first choosing the antinodal direction to be orthogonal to the plane of the interface (|] , second by restricting 
electron transfer along the normal of this interface. Nevertheless, the confrontation of our theoretical results with 
transport experiments in high-T c compounds may prove too difficult, as surface defects are likely to be present in 
such compounds. Such defects are responsible for a reduction of the order parameter in the vicinity of the junction. 
Moreover, the presence of Andreev bound states - also induced by the impurities at the interface - renders the present 
approach inadequate. Furthermore, surface roughness provokes a randomization of the tunneling directions so that 
the above scenario for antinodal tunneling loses its practicality. 

Previous theoretical works addressing the role of phase fluctuations (24j| on Andreev reflexion in the pseudogap 
phase focusscd on the mildly underdoped case: the hole density is chosen on the right hand side of the zero 
temperature insulator-superconductor quantum critical point (QCP). In this regime quantum fluctuations may renor- 
malized the superfluid density but the physics turns out to be essentially that of the two-dimensional classical XY 
model (renormalized classical regime in the spirit of the n = 2, non- linear cr-model). Transverse superfluid velocity 
fluctuations due to the motion of vortices are expected to dominate. Because the time scales for the vortex motion, 
as measured by THz spectroscopy [ p5| are much larger than the Andreev time ta = ^/A, the correlation of the phase 
was treated as static (on the scale of ta)- The phase correlator decays in space over a lengthscale £ c / which is the 
classical correlation length Jl6| . If tunneling occurs on an area which is much smaller than £ c ; , then an Andreev signal 
is expected. 

Our calculation (albeit in 1 dimension), focuses on quantum fluctuations. Fluctuations are included in a non- 
perturbative and systematic manner, regardless of the phase propagator which is used in the end. As a general trend, 
quantum fluctuations were found to have a tendency to reduce the Andreev reflexion signal. Note that other models 
of pseudogap also rely partly on quantum fluctuations p(| , while others incorporate aspects of ID physics (stripes) 
. Extensions of our calculation beyond ID - together with surface roughness - or with competing orders could be 
envisioned to achieve a more realistic description of STM experiments on high-T c compounds. 
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